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Abstract 

Precise /3-branching-ratio measurements are required in order to determine /t-values as part of our program to test 
the Electroweak Standard Model via unitarity of the Cabibbo-Kobayashi-Moskawa matrix. For the measurements 
to be useful in this test, their precision must be close to 0.1%. In a branching-ratio measurement, we position the 
radioactive sample between a thin plastic scintillator used to detect /3-particles, and a HPGe detector for 7-rays. 
Both (3 singles and /3-7 coincidences are recorded. Although the branching ratio depends most strongly on the HPGe 
detector efficiency, it has some sensitivity to the energy dependence of the /3-detector efficiency. We report here on 
a study of our /3-detector response function, which used Monte Carlo calculations performed by the Geant4 toolkit. 
Results of the simulations are compared to measured /3-spectra from several standard /3-sources. 

Key words: Precise /3-branching-ratio measurements; the Electroweak Standard Model; unitarity of the 
Cabibbo-Kobayashi-Moskawa matrix; Monte Carlo simulations; Geant4 efficiency calculation for plastic /3-detector. 



1. Introduction 

Knowledge of the total efficiency of a plastic /3- 
detector is crucial to our precision experiments test- 
ing the Electroweak Standard Model. We measure 
the ft values for superallowed 0+ nuclear 
transitions, from which we obtain the value of Yudi 
the up-down quark-mixing element of the Cabibbo- 
Kobayashi-Moskawa (CKM) matrix. This requires 
that half-lives, branching ratios and decay energies 
all be measured with high precision, 0.1% or bet- 
ter. (The most recent complete review of this work 
can be found in Ref. [1], with an update in Ref. 
[2]). Since branching ratios are typically determined 
from the intensities of /3-dclayed 7 rays that are ob- 
served in /3-7 coincidence measurements, it is the 7- 
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ray detector's efficiency that is most crucial in aim- 
ing for 0.1% precision. However, a good knowledge of 
the energy-dependence of the /3-detector's efficiency 
is also required. Here we report studies of our /3- 
detector's response function, with source measure- 
ments and Monte Carlo calculations performed with 
the Geant4 (version 4.9.0) toolkit [3]. 

In a typical measurement of a /3-decay branching 
ratio (see, for example, Ref. [4,5]), we implant a ra- 
dioactive species into Mylar tape, then rapidly move 
the tape to a shielded counting station, where the 
sample is positioned between a 1-mm-thick plastic 
scintillator to detect /3-particles, and a 70% HPGe 
detector to detect 7-rays. We record both (3 singles 
and (3-"f coincident events. Since the efficiency of 
the HPGe detector has been very precisely deter- 
mined [6] , to first order a /3-branching ratio is simply 
given by the measured number of coincident 7 rays 
that follow that (3 branch divided by the total num- 
ber of /3 singles from all branches; the /3-detector 
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Fig. 1. Schematic drawing of the /3-detector. Labeled dimen- 
sions are given in Table 1; others appear in the text. Inci- 
dent /3-particles from a radioactive source pass through a 
thin havar-foil window (3) and into the plastic-scintillator 
disc (2). The disc is recessed into a Incite light guide (4), 
which is coupled to a phototube (not shown in the figure) 
at surface 5. The scintillator and light guide are enclosed in 
a cylindrical cover (1) made from PVC. 

efficiency simply cancels out and need not be deter- 
mined. For high precision, though, it becomes nec- 
essary to account for the slightly different efficiency 
of the P detector for each f3 transition, a difference 
that affects the measured intensities of the coinci- 
dent 7 rays. 

The energy dependence of our /3-detection effi- 
ciency is caused principally by the low-energy elec- 
tronic threshold, which removes a slightly different 
fraction of the total (3 spectrum for different end- 
point energies. Since our threshold is at ~80 keV 
and, for our superallowed-decay studies, end-point 
energies are typically 2 MeV or more, we lose at 
most a few percent of the total /3 particles for any 
single transition. Thus the change in this loss from 
transition to transition in the same decay is even 
smaller. Nevertheless, the precision we strive for in 
our branching-ratio measurements is very high and 
we seek to account reliably for that energy depen- 
dence. Even though a Monte Carlo code like Geant4 
should be well suited to simulating the response 
function of a thin plastic scintillator, we consid- 
ered it important to test and evaluate the code's re- 
sults first, by comparing them against experimen- 
tal data taken with several f3 sources that also emit 



conversion-electrons: ^^■^Ba, ^■^^Cs and ^°^Bi. 

The Geant4 toolkit used in treating the trans- 
portation of P and 7 particles through matter is 
both modular and flexible, especially in the descrip- 
tion of low-energy electromagnetic processes down 
to 250 eV (see, for example Ref. [7,8]). In addi- 
tion, it is also possible to simulate rather compli- 
cated 3-D geometry, select a variety of materials 
and decay products (including radioactive ions), and 
choose how to handle the physical processes govern- 
ing particle interactions. Moreover, it provides out- 
put of the simulated data at different stages in the 
calculation and under various selection criteria. An 
overview of recent developments in diverse areas of 
this toolkit is presented in Ref. [9]. 

2. Detector Arrangement and 
Measurements 

The /3-detector assembly is illustrated in Fig. 1, 
with detailed dimensions given in Table 1 . It consists 
of a 1-mm-thick Bicron BC404 scintillator disc re- 
cessed into a cylindrical Lucite light guide, to which 
it is optically coupled. The light guide, in turn, is op- 
tically coupled to a photomultiplier tube (not shown 
in the figure and not included in the simulations). 
Optical cement (BC-600) from Bicron was applied 
to the surfaces between the scintillator and the light 
guide, and between the light guide and the pho- 
tomultiplier tube (R329P from Hamamatsu). The 
scintillator disc, light guide and the last 13 mm of 
the phototube are enclosed in an opaque cylindri- 
cal shield made from 1.5-mm-thick polyvinyl chlo- 
ride (PVC). The opening in the scintillator end of 

Table 1 

Measured detector dimensions used in our Monte Carlo cal- 
culations. Letters in the second column correspond to labels 
in Fig. 1. 



Detector parameter 




Value (mm) 


Outer shoulder, 


Di 


63.50 


Inner shoulder, 


D2 


60.33 


Light guide, 


D3 


49.20 


Outer PVC cover, 


D4 


55.63 


Hole in PVC cover, 


Ds 


38.61 


Plastic scintillator, 


De 


38.10 


Length of PVC cover 


Li 


53.85 


Length of outer shoulder 


La 


39.62 


Length of light guide 


L3 


39.24 


Length of inner shoulder 


U 


12.70 
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the PVC shield is sHghtly larger in diameter than 
the scintillator disc, and is covered with a pin-hole- 
free, 5-/xm-thick havar foil. The /3-particles enter the 
detector assembly through this foil with essentially 
negligible energy loss. 

For our measurements of the detector response 
function, each radioactive source was placed at a 
distance of 13 mm from the havar window of the 
detector assembly and was axially aligned with it. 
The distance was determined with the aid of an Ac- 
cuRange 600'^'^ Laser Displacement Sensor (model 
AR600-4000) [10], which measures distance with an 
absolute precision better than 0.1 mm. Both detec- 
tor and source were placed on stands on a lab ta- 
ble (in air at atmospheric pressure) as far away as 
possible from other objects. We used three differ- 
ent 37-kBq sources - ^^^Ba,, "^Cs and ^^^Bi. All 
three were open sources sold by Isotope Products 
Laboratories as "conversion electron sources." Each 
source, being specially prepared to minimize scat- 
tering or degradation of the emitted electrons, had 
been deposited as a 5-mm-diameter spot on a thin 
foil - stainless steel with a thickness of 51 fim in the 
case of ^"^Bi, aluminized Mylar with a thickness of 
6 /im for the other two - and covered only by a 100- 
fig/cm^ acrylic film. The source-holder geometry is 
shown in Fig. 2. 

All three of these radioactive sources emit 7 rays 
as well as electrons and, although our thin scin- 
tillator is relatively insensitive to the former, we 
nonetheless took extra precautions to ensure that 
we were only studying the detector's response to the 
latter. In addition to recording a spectrum from each 
source as already described, we also recorded a sec- 
ond spectrum with a 2-mm-thick aluminum plate 
inserted between the source and the detector. This 
plate was thick enough to remove all the f3 par- 
ticles without significantly attenuating the 7 rays. 
We then subtracted this second spectrum from the 
first, and considered the resultant spectrum to be 
a "pure" /3 spectrum. This method has one flaw, 
however: the spectrum obtained with the aluminum 
plate includes some contribution from the bremm- 
strahlung created by the /3 particles as they stop. 
Thus, we took the same approach with the calcu- 
lated Monte Carlo spectra: we made two calcula- 
tions for each source, one with an aluminum plate 
and one without, took the difference between them 
and then compared that difference spectrum with 
the "pure" experimental /? spectrum. 
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Fig. 2. Radioactive-source geometry: the support ring is 
made from aluminum. 

3. Geant4 Physics Model 

The Geant4 Simulation Toolkit includes a series 
of packages for the simulation of electromagnetic in- 
teractions of particles with matter, specialized for 
different particle types, energy range and specific 
physics model. In our work, we considered only elec- 
trons, 7 rays and x rays, and used three different 
physics models for the electromagnetic (EM) pro- 
cesses: the standard EM package, the low-energy 
EM package and the Penelope EM package. In all 
cases, fluorescence emission, Rayleigh scattering and 
Auger interactions were included in the EM physics 
model where appropriate. The standard EM package 
is based on an analytical approach [11,12,13]; its ef- 
fective energy range is nominally between 1 keV and 
100 TeV but it neglects atomic effects and is mainly 
used in high-energy physics applications. The low- 
energy package is optimized for our energy region 
and extends the range of validity for electrons and 
photons down to 250 eV [7,8] and even below. The 
Penelope package is an alternative low-energy im- 
plementation; it is a re-engineered version [14] of 
the original PENELOPE Monte Carlo code [15,16]. 
Detailed information on these packages, as well as 
on the design of the Geant4 toolkit can be found 
in Ref. [3] and in the Physical Reference Manual of 
Geant4 [17] and references therein. 

The validation of the Geant4 electromagnetic 
physical processes is important in order to reach 
an adequate level of precision in applications such 
as ours. Systematic and extensive validation is an 
on-going process in the Geant4 collaboration and 
recently microscopic quantities such as cross sec- 
tions, angular/energy distributions, attenuation 
coefficients, stopping powers and ranges have also 
been examined in a systematic way [18], and their 
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Fig. 3. Upper panel: Energy deposition in the /3-detector as generated by Geant4 Monte Carlo calculations for the decay of 
^^^Cs. Three different EM physics models are used and all three are plotted although they are so similar that they cannot be 
distinguished. Lower panel: To clarify the differences among the results for the three EM physics models, the low-energy EM 
package is used as our standard, and differences between each of the other two models and the low-energy package are plotted as 
residuals in standard-deviation units. The dashed line corresponds to the Penelope EM package (compared with the low-energy 
package); while the dotted line corresponds to the standard EM package (also compared with the low-energy package). 



compatibility with reference data from the U.S. 
National Institute of Standards and Technologies 
(NIST) has been estabhshed [14]. 

Here we have restricted ourselves to a compar- 
ison between Geant4 calculations and experiment 
for the energy deposited by electrons in a thin plas- 
tic scintillator, where we have used the very simple 
laboratory geometry already described so that the 
Monte Carlo geometry could reproduce it exactly. 
As can be seen in the top panel of Figure 3, under 
these conditions the three physics models - stan- 
dard EM, low-energy EM and Penelope EM - gen- 
erate energy spectra that differ very little from one 
to another. The bottom panel of the figure shows 
the normalized residuals between the first and sec- 
ond models, and between the second and third ones: 
there are small but perceptable differences below 200 
keV but nothing significant above that energy. We 
also compared the total /3-efficiencies obtained from 
the three EM physics models. Including all ener- 
gies between 50 and 1000 keV, the calculated Monte 
Carlo efficiencies were 15.16(3)%, 15.18(3)%, and 
15.11(3)%, respectively. If the low-energy thresh- 



hold was increased to 75 keV the calculated efficien- 
cies were 13.79(3)%, 13.84(3)% and 13.82(3)%. For 
both thresholds the three physics models yielded 
statistically identical results. In short, for our pur- 
poses we find nothing in these results to choose be- 
tween the three available physics models. 

Even so, we chose to use the low-energy EM model 
in all the Monte Carlo simulations presented in the 
remaineder of this paper. Although it took consid- 
erably more computer time per calculation than did 
the standard EM model, we considered that it was, 
in principle, more appropriate to our energy region 
since it was specifically designed for better perfor- 
mance at low energies. A detailed inter-comparison 
of results from the three physics models for the total 
efficiency of the plastic /?-detector as a function of /? 
energy will appear in a subsequent publication [19]. 

4. Geant4 Geometry 

In defining the laboratory geometry in Geant4 we 
included all the components of the detector assem- 
bly and source housing (see Figs. 1 and 2, and Ta- 
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Table 2 

Composition of the different materials used in the Monte Carlo simulations performed in this work. The tabulated values 
correspond to the element mass fraction in each material, given in percentages. Material densities are also given. 



Chemical 
element 


Air 


Acrylic 


Mylar 


Havar 


Stainless 
Steel 


PVC 


Lucitc 


BC404 


Aluminum 


H 




0.71 


4.20 






4.84 


8.07 


8.45 




C 


0.01 


8.52 


62.50 


0.04 




38.44 


59.97 


91.55 




N 


75.53 


90.77 
















O 


23.18 




33.30 








31.96 






Ar 


1.28 


















Be 








0.01 












CI 












56.73 








Si 










1.00 










Cr 








17.04 


19.00 










Ni 








12.50 


10.00 










Fe 








16.34 


68.00 










Co 








41.04 












Mo 








3.14 












Mn 








1.45 


2.00 










Al 


















100.00 


W 








8.44 












Density 
g/cm^ 


0.0012 


1.190 


1.390 


8.300 


8.020 


1.380 


1.185 


1.032 


2.700 



ble 1), with everything placed in air. Special care 
was taken to include all components of the various 
materials, with the natural isotopic abundances for 
each element properly accounted for. Table 2 lists 
the composition and density of all materials used in 
the Geant4 calculations. 

As explained in Section 2, two Monte Carlo sim- 
ulations were performed for each source, one with a 
2-mm-thick plate of aluminum located between the 
source and the detector, and one without it. The 
spectrum we plot for each source is the difference be- 
tween these two spectra. In defining the geometry in 
Geant4 we incorporated a plate and chose its mate- 
rial to be either aluminum or air, depending on the 
desired effect. Fig. 4 shows the tracks of electrons 
for both types of simulations. 

5. Geant4 Parameter Control 

Geant4 allows the user to define different regions 
in the experimental setup, and to set a different 
particle-production threshold in each one [9]. This 
capability allows for simulation accuracy and speed 
optimization according to the needs of a particular 
experiment. We defined three regions: the radioac- 
tive source, the thin plastic scintillator and every- 
thing else. In the first two, the threshold for produc- 



ing secondary particles was kept very low in order 
to simulate all physics interactions as well as possi- 
ble; in the third, a much higher production threshold 
was chosen. This approach considerably reduced the 
computing time without compromising the results. 
A test run with low thresholds in all regions did not 
reveal any significant differences from the sped-up 
version with different thresholds in different regions. 

It was pointed out by Kraev [20] that, to get 
good agreement with experiment, it is important to 
choose two parameters particularly carefully in the 
low-energy EM physics model of Geant4: the cut for 
secondaries (CPS) value, which determines the pro- 
duction threshold for secondary particles, and the 
/^-parameter, which limits the step size for tracking 
/3-particles at a material boundary. In this work, we 
used CPS — 10 /im, which is recommended as an ac- 
ceptable compromise between a "good" description 
of the scattering processes and a reasonable com- 
putation time. For the /^-parameter we used 0.02 
as recommended in [20]; this is actually the default 
value for the Geant4 version 4.9.0. 

6. Comparison with Experiment 

To simulate the decay of a radioactive nuclide 
with Geant4, it is possible to define each 7 transi- 
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Fig. 4. Schematic view of the geometry arrangements for 
Monte Carlo simulations, showing sample electron trajecto- 
ries. The upper panel shows trajectories when there is only 
air between the source holder (right) and the detector hous- 
ing (left). The lower panel shows them when a 2-mm-thick 
aluminum plate is introduced. 

tion, internal-conversion line and /?-decay spectrum 
individually and require Geant4 to transport all 
particles through the specified materials and de- 
termine the spectrum in the scintillator. However, 
the code also ofi^ers a radioactive-decay module, 
which generates all the decay components radiated 
from a specified source using information extracted 
from the Evaluated Nuclear Structure Data File 
(ENSDF) [21]. 

To generate the Monte Carlo emission spec- 
tra we began by programming Geant4 based on 
the radioactive-decay module. The primary elec- 
tron spectrum emitted from ^"^^Bi generated with 
the radioactive-decay module is shown in Fig. 5. 
When repeating this procedure for ^'^'^Ba, to our 
surprise we found that the electron emission spec- 



trum produced by the radioactive-decay module 
of Geant4 was simply not correct, yielding relative 
conversion-electron intensities in significant dis- 
agreement with ENSDF data. The emission spec- 
trum from ^^^Cs also turned out to be incorrect, 
but here the main problem was more subtle: there 
are two /3-decay branches from ^■^''Cs, which are 
both treated by Geant4 as allowed. In fact both 
transitions are forbidden, with shape-correction 
factors that have been determined by Behrens and 
Christmas [22] from experimental data. In addition, 
the radioactive-decay module gives the incorrect 
intensity for one of the conversion electron lines of 
"7Cs (/655.7 kov=l-10% instead of 1.39%). In both 
these decays - of ^^'^Ba and ^'^''Cs - we bypassed 
the radioactive-decay module and inserted each 
decay mode and transition individually, with the 
correct intensities for the conversion electrons and 
the correct shape for the forbidden /3 transitions. 
For these we used the General Particle Source mod- 
ule available in the Geant4 [23], which allows the 
user to define standard energy, angle and space 
distributions of the primary particle. 

Based on a primary spectrum thus generated for 
each source, the Monte Carlo code then determined 
the total energy deposited in the scintillator. How- 
ever, before this result could be compared with the 
experimental spectrum, it was necessary to add 
the effects of statistical fiuctuations introduced by 
the processes of light production and transmission, 
as well as photomultiplication and electronic pulse 
analysis. For this purpose, we looked to a published 
study of the response of a plastic scintillator to 
mono-energetic beams of positrons and electrons 
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Energy, keV 

Fig. 5. Decay spectrum for ^'^'^Bi generated by Geant4 with 
its internal radioactive-decay module activated. Only elec- 
trons are shown. 
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Fig. 6. In the top panel, the measured spectrum (thick 
sohd Une) for the decay of ^^^Ba is compared with the 
Geant4-simulated result (thin dashed hne). The low-energy 
EM package was used. Residuals in standard-deviation units 
are plotted in the lower panel. The reduced in the energy 
range 80 - 406 keV is 0.4. 

[24], which tabulated the width of the full-energy 
Gaussian peak as a function of energy between 0.8 
and 3.8 MeV. Since we also needed to deal with 
energies lower than that, we took the width to be 
linearly dependent on energies below 0.8 MeV. 

Our procedure was to take the scintillator spec- 
trum produced by Geant4 and process it by a ran- 
domization algorithm written in C-f -f in the ROOT 
[25] analysis framework. In essence, this process 
spread the number of counts in each energy bin into 
a Gaussian distribution centered at the original 
energy and with a width, a, taken or extrapolated 
from Ref. [24]. The results could then be compared 
directly with the measured spectra. 

The measured spectra for all three sources were 
taken under exactly the same conditions. A descrip- 
tion of the data acquisition system used in these 
measurements is given in ref. [4] and references 
therein. The amplifier gain, photomultiplier high 
voltage, and low-energy electronic threshold re- 
mained unchanged for all three sources. The thresh- 
old was chosen to be identical to the setting used 
during our branching-ratio measurements of super- 
allowed Fermi /3-decays. (This hardware threshhold 
was higher than the software threshhold.) This en- 
sured that the results of our comparisons between 
source measurements and Monte Carlo simulations 
could be readily apphed to our accelerator-based 
measurements. 

In comparing our measured spectra with Monte 
Carlo results, we the slope of the energy calibra- 
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panel, the measured spectrum (thick 



solid line) for the decay of ^^'^Cs is compared with the 
Geant4-simulated result (thin dashed line). The low-energy 
EM package was used. Residuals in standard-deviation units 
are plotted in the lower panel. The reduced in the energy 
range 80 - 697 keV is 4.0. 

tion (i.e. the energy per channel) and its offset ((i.e. 
the zero-energy channel) corresponding to each mea- 
sured spectrum as adjustable fit parameters, which 
were used to optimize the agreement with the Monte 
Carlo spectrum. We wrote our own C-|— I- ROOT 
program to accomplish this purpose. The resulting 
comparisons for our three sources, as well as the nor- 
malized residuals, for ^'^^Ba, ^"^^Cs and ^"^Bi, appear 
in Figures 6, 7 and 8 respectively. 

The agreement between the Geant4 Monte Carlo 
simulations and experiment is good for all three 
sources considered in this work, ^■^^Ba, ^"^^Cs and 
^°^Bi, although for the latter one the normalized 
of 8.9 is less impressive than for other two. In that 
case, the energy range, which extended from 80 to 
1143 keV, was much greater than for the other too. 
This may be partly responsible for the higher 
but another possibility is that our simple linear ex- 
trapolation of the results in Ref. [24] does not fully 
describe our system's response function at low en- 
ergies. If in fact the peak resolution were somewhat 
worse than this extrapolation indicates - a not un- 
reasonable possibility ~ then the agreement with ex- 
periment would be considerably improved. It is also 
worth noting that the slight non-allignment in peak 
positions, which is evident in all three spectra, can 
be explained by possible small non-linearities in the 
experimental energy response. 

Since the amplifier gain and photomultiplier high 
voltage were kept the same for all three measure- 
ments, the fitted slopes and offsets obtained for all 
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Fig. 8. In the top panel, the measured spectrum (thick 
sohd line) for the decay of -^"^Bi is compared with the 
Geant4-simulated result (thin dashed line). The low- energy 
EM package was used. Residuals in standard-deviation units 
are plotted in the lower panel. The reduced in the energy 
range 80 - 1143 keV is 8.9. 

three spectra should have been very nearly the same. 
In fact they were, but, as a quantitative measure of 
consistency, we took the average values for the slope 
and offset, and again evaluated the normalized 
for the comparison between experiment and simu- 
lation. The new values for the normalized were, 
of course, somewhat increased, being 3.6, 9.5, and 
16.3 for i33Ba, ^o^Bi, and ^^^Cs, respectively, but the 
agreement is still quite satisfactory. Most crucially 
for our purposes, the low-energy thresholds in all 
three cases were in close agreement, with 80±3 keV 
being the common value. 

As stated in the introduction, the result of the 
simulation that matters most to us is how well it 
reproduces the fraction of the total P spectrum 
that lies above some low-energy threshold, typi- 
cally around 80 keV. The spectral details at higher 
energy are only important to us to the extent that 
they change the fraction of counts recorded above 
threshold. The fact that the fits performed with 
the low-energy threshold as a free parameter give 
the same result for all three isotopes, serves as im- 
portant assurance that the fitted threshold value is 
consistent with the actual one. 

Although the activity of the radioactive sources 
that we used for this work are nominally 1 /iCi (37- 
kBq), the accuracy of this value was only quoted 
to and approximate ±15% by the supplier. So that 
we could get a more precise value for our /3-detector 
efficiency, we made our own measurement of the 
^"^^Bi source activity using a well-calibrated HPGe 



7-detector [27] to detect the known 7 rays from the 
decay. In this way we established the source activity 
to be 1.31(1) /iCi at 16 January, 2008. Now, knowing 
the activity of the source as well as the low-energy 
detection threshold already obtained from our fit, 
we could deduce from our experimental data the ab- 
solute efficiency of the /^-detector to be 3.48(2)% at 
the distance of 13.2(1) mm. With exactly this ge- 
ometry, the Geant4 simulation yielded an absolute 
efficiency of 3.50(1)%, in excellent agreement with 
experiment . 

7. Conclusion 

The electron spectra we have obtained from 
Monte Carlo simulations are generally in good 
agreement with experimental data. All identifiable 
features that are present in the experimental spec- 
tra are reproduced in the simulated ones, and in 
most cases their relative intensities agree as well. 
Furthermore, the threshold energy and absolute 
efficiency are well reproduced. Our results clearly 
demonstrate that Geant4 version 4.9.0 can be used 
effectively to simulate the /3-spectra as measured 
by our thin plastic scintillator detector. 

In particular for our application, where we need 
only rely on the simulation to determine the energy 
dependence of the /3-detector 's total efficiency with 
a low-energy threshold at ^^80 keV, it is clear that 
Geant4 will provide the precision we require. 

However, we have also shown that the radioactive- 
decay module included in Geant4 to simulate the ini- 
tial radiation from a radioactive source should only 
be used with care: it is essential to check its output 
carefully. In particular, the intensities of conversion- 
electron lines produced by this module were some- 
times found to be incorrect and, if any forbidden (3- 
decay branches are involved, their spectrum shapes 
may not be correctly generated. 
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